function res = Romberg(fun,xl,xr,n,N)
% ROMBERG - Romberg 方法求积分
%   
if N ==1,
    res = multitra(fun,xl,xr,n);
    return
end

res = (Romberg(fun,xl,xr,2*n,N-1)-4^(1-N)*Romberg(fun,xl,xr,n,N-1))/(1-4^(1-N));
    
end
